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In strongly magnetized backgrounds, electrically charged vector bosons can condense. When 
Qh considering the p-meson sector of QCD, this results in the existence of a superconducting state 

with Abrikosov vortices. As a next order effect, also the neutral mesons condense and one expects 
the resulting lattice to be a mix of neutral and charged condensates, opening a possibility to the 
existence of a superfluid property. We will show how the detailed structure of the vacuum can be 
calculated and which implications it has on the nature of the exotic superconducting state. 
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Electroweak versus QCD-vacuum in strong magnetic fields 



1. Introduction 

It was proposed last year [1] that in a strong magnetic field the vacuum may become an electro- 
magnetic superconductor. For very strong magnetic fields, the lightest of the charged vector mesons 
are likely to condense in a similar way as the electroweak vector bosons [2]. The p-mesons have 
a polarization with the right quantum numbers to see its energy reduced by an increasing magnetic 
field, a property that is not shared with the lighter 7T-mesons. Hence the validity of the restric- 
tion to the p -meson sector and the stability of the condensates, as ^-mesons are no longer favored 
energetically. 

While it is a hope to observe the effects of such condensates experimentally, the required 
static magnetic field is way beyond experimental capabilities. However, a beam collision with 
very small impact parameter could create a magnetic field with the required strength [3], about 
10 16 r. This setting is promising, but does share the absence of electrical fields with the theoretical 
approximation. It is unclear whether the new vacuum could still appear in this setup and if it 
could create a detectable signal. If not, the relevance may be in Big Bang scenarios, as the early 
universe may have been strongly magnetized [4] and have suffered some effects of the exotic state 
this implies, as we will show. From a theoretical point of view there's also a strong incentive 
to study this effect, as the results seem to have a broader application base than the p-mesons 
only. For example, W-condensation in the electroweak model [2] can be better understood by the 
very same -or at least similar- calculations and in particular one finds that the condensates imply 
superconductivity [5]. The obtained solution to the equations of motion joins other well known, 
non-trivial such solutions like monopoles, cosmic strings and sphalerons, all interesting in their 
own right[6]. 

Aside from obtaining a better picture of the vacuum configuration and its superconducting 
nature, the calculations presented give some insight in the potential superfluid structure of the vac- 
uum. While no proof of superfluidity was found, and in particular the superfluid London equations 
are not satisfied, we find a substructure in the Abrikosov lattice consisting of 'superfluid' vortices 
[7]. These vortices give a non-trivial winding of the phase of neutral fields, just like in superfluids, 
but it remains unclear whether our exotic medium is a superconducting/superfluid hybrid or just a 
very peculiar superconductor. The superconductivity, at least, was confirmed in numerical lattice 
calculations [8] and also hinted from holographic approaches [9]. 

As a last interesting fact, we found the superconductivity to be strongly influenced by the 
neutral condensate such that, in the core of the Abrikosov-vortices the London-equation is still 
satisfied, unlike the ordinary case where it vanishes at the core. This means in some sense the 
superconductivity is stronger than in the ordinary case. 

In these proceedings we will go into detail on how the Abrikosov lattice can be obtained 
starting from nothing but the Lagrangian of the theory. It corresponds roughly to the paper [7], 
work done in collaboration with Maxim Chernodub and Henri Verschelde. 

2. Effective potentials 

To find an effective potential that can give a description of the Abrikosov lattice we will ba- 
sically need three steps: we describe the model in the transversal (to the external magnetic field) 
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plane, then make a small field approximation and eventually solve the simplified equations of mo- 
tions and write the energy density as a function of one scalar field only. We will repeat the steps 
for three interesting cases: the ordinary Ginzburg -Landau model for comparison, the electroweak 
theory studied by Ambjorn and Olesen and the p-meson sector of the QCD-vacuum, described by 
the DSGS model [10]. Many of the calculations can be found already in [1] and more recent [7]. 

2.1 Ginzburg-Landau model 

The Lagrangian we start from is given by 

L=\D^\ 2 -V(\$\ 2 )-^F 2 V ; (2.1) 

V{x)= l -X{x-$f; D ti <j> = d fl (j>-ieA fl <j>. (2.2) 

With the identification 

= e x +i# 2 ; d = G\- iff 2 ; (2.3) 

for any operator, we can rewrite this as an energy density, putting all longitudinal ((0, 3)-directional) 
dependence to zero and lowering all (spacial) indices. 

£=\D$\ 2 + eF i2 \$\ 2 + V(\4)\ 2 )+ l -F? 2 . (2.4) 

Here surface terms are neglected and we used [D,D]<p = —2eFi2<j>- 

Now we are ready to make an approximation near the critical point for a sufficiently high 
external field eB = e(Fn) = j^>q, where all terms in the above expression become positive. 

eF 12 >^<l>$^<j>=0 (2.5) 
eFn < \§l => 50 = G^ 2 ) ~ (2.6) 

Neglecting the contribution of the covariant derivative in a well-chosen gauge as well as a period- 
icity requirement for the lattice, gives an Abrikosov ansatz [11] for test-space of solutions for 0, 
that will be used for a variational method: 

D0„ = Hr - e -B^ n = =► *, = e-iWe-to't-i'*-*, (2.7) 

with Lb = y|§. It is easily seen that with this form translations in real and imaginary directions 
give (n,N G N): 

(j) n (z + iL B ,z- iL B ) = e r ^ <j) n (z,z) (2.8) 
Uz + NL B ,z + NL B ) = e iK '*B (j> {n+N) (z,z) (2.9) 

Neglecting the n independent gauge transformation, we find a periodic solution 

<Kz,z)=£c„<Mz,f) (2-10) 
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if C n+ N = C n for some N and for all n. Fixing N, we now have an ,/V-parameter ansatz to be plugged 
into the energy functional. 

Lastly we determine the appropriate energy from (2.4). The first term can already be dropped, 
as our test function makes it negligible. The equation of motion for the gauge field gives further- 
more: 

d(F n + e\<j)\ 2 ) =0 (2.11) 
which implies the argument is constant and equal to its average value, thus: 

F n =B + e(\f\)-e\$\ 2 . (2.12) 

The last step is straightforward: 

(,f) = (eS-^ o 2 )(|0| 2 ) + ^(|^| 2 ) 2 + (^-^)(|0| 4 ) + o 4 (2.13) 

Plugging in (2.10) allows a numerical minimization in the iV-dimensional parameter space. 

Interestingly, the third term in the above expression can be rewritten in the vacuum solution 

as: 

/M^-yXm 2 ) 2 (2-14) 

with jSa the Abrikosov-parameter determining the lattice packing, unlike (|0| 2 ) that is independent 
of it. When = this ordering becomes degenerate, indicating a phase transition. Indeed, this 
corresponds to the Bogomolnyi limit at which type II superconductivity goes to type I and vice 
versa. 

A last point we can look at is an explicit proof that we have superconductivity. Therefore we 
need the covariant form of the London-relation 

df n s e 2 : 

— = -^E' (2.15) 
at m 

with E the electrical field and a the conductivity, essentially depending on the density of cooper- 
pairs. We can determinge equations of motion assuming a very small electrical field in the longi- 
tudinal direction (parallel to the external magnetic field, but all other components still transversal. 
We get for the current = d v F vll : 

d J 3 -d 3 J = 2e 2 \<p\ 2 E 3 (2.16) 



2.2 Electroweak vacuum 

Now that we have gone through the simplest case in much detail, we can apply this to a more 
complicated model. The Lagrangian of the electroweak model in the unitary gauge reads: 

L = - l -(W« v ) 2 -^X^ + lD^l 2 -V(\(t>\ 2 ), D^ = d^-i^£^Z^, (2.17) 

with the photon field = cosdX^ + sin0W^, Z M = -sindX^ +cos0W^ and W~ = (W* + 
iW 2 
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Taking the convention W = W l + iW 2 , neglecting the component W and gauging to the unitary 
gauge, the expression for the energy density for transversal components only is: 

^ = \ F h + \\^W\ 2 + ^( 8 -^-eF i2 -gcosdZ l2 )\W\ 2 + 8 j\W\ 4 

+ ^ 2 + 4Se |z|2 + l ^ |2 + y(l</,|2) (2 " 18) 

with = d W + ieAW + ig cos 6ZW. Knowing that in the ordinary vacuum gives a mass to the 
W, we observe the same compensating effect of the magnetic field in the mass term as we did in 
the ordinary superconductor. Very important is however that now the sign of this term is reversed, 
a direct consequence of the spin-1 nature of the IV-field. It is easy to calculate that for positive F\2 
this effect is absent for the W component and that it will not contribute a vacuum expectation value. 
Leaving it out was thus justified. 
Now we have 

eF 12 <y0 o 2 ^W = O (2.19) 

eFn> y0o 2 ^@W = 0{W 2 )~O (2.20) 

We can in fact assume all terms to be the same small order &(W 4 ), as they are in an approximation 
up to &(W 2 ) all still vanishing. With Z ~ <p — <po ~ &(W 2 ), we can solve their respective equations 
of motion and the one for the photon field and find: 

z =i cos6d jdh^J w ^ ^°=^^4 |M/|2; ^=*^\ 2 ->\ 2 y> 

, 2 ' ' (2-21) 

with M 2 = and M 2 =X^) 2 . This is very similar to the Abelian case, but with a different sign 

for \W\ 2 , indicating the magnetic field and the superconducting condensate are in fact proportional 
to each-other, in stead of the usual inverse correlation. This property is called anti-ferromagnetism 
and it indicates that the superconducting state, rather than being destroyed at a sufficiently high 
magnetic field, is created at a certain critical field strength. 

By assuming the very same ansatz (2. 10) for W as before, we can find a variational lattice 
solution if we rewrite the energy functional in terms of the W-condensate only. This involves some 
careful bookkeeping in coefficients, but is fairly straightforward: 

<*> = £<|*ry + 1 ( a4 -es^wf) + IbL + d*« ^ _ _J^) |«f 

(2.22) 

Clearly, the Bogomolnyi limit appears for Mz = Mh, exactly as was already shown in Ambjorn and 
Olesen's original work. 

2.3 QCD vacuum 

Armed with insight from the previous models, we can put the machinery to work in case of an 
effective model for the p-meson sector of QCD. As argued in the introduction, this will hold the 
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dominating contribution to a vacuum condensate in a strong magnetic field. The DSGS Lagrangian 
is given by [10]: 

L = -\{PI V ) 2 ~ \f^ v + ^(p«) 2 + ^p£V v - (2-23) 

For consistency, we will assume the same notation for the charged mesons as we did previously 
for the vector bosons when rewriting the theory in the transverse plane. It is good to know that the 
convention in the original paper [1] is found by p i— > 2p in the formula below. 

Neglecting the p -polarization in favor of the p for the same reason as we dropped W previ- 
ously, but the other way around due to the conventions in [10, 1]. We obtain for the transversal 
degrees of freedom: 

<? = \Fh + \m\ 2 + \ (4 +ft /g> - eFn) |pf + f |P| 4 + \{f^f + i*><°>| 2 - £iWg> 

(2.24) 

with @p = dp -ieAp + ig s p^ 0) p and = -|(<9p (0) - dp (0) ). The similarity with the elec- 
troweak model is obvious and we proceed in exactly the same way, making the appropriate estimate 
for all fields assuming e(F\2) = eB> m\. The equations of motion give us 



P" 

P (0) = ~i 8 i d^- 2 \p\ 2 ; F l2 =B+ e -\pf -U\p\ 2 )- V^J^-M' ^ 
^do—nig I 2. I da—iriQ 

using (/I? ) = and n& = m 2 D /{\ — ^). This leads tediously to 

{*) = l&ext + j (IPl 2 ) 2 + \{m 2 p -eB ext ){\9\ 2 ) + f rn 2 p / \P\ 2 _ d ^ +n ^ \P\ 2 ) » (2-26) 

where we replace p by the Abrikosov ansatz illustrated before and proceed to the numerical mini- 
mization. 

It must be noted that the coefficient of the last term is much larger than in the electroweak 
case, as there's no compensating second term in it. This makes the algoritm much more stable, and 
is exactly the reason we obtained so far only results for this case, the electroweak structure should 
soon follow [5], 



3. Results 

The results show an ordering of vortices in a hexagonal structure, which can be obtained in 
a 2-parameter Abrikosov ansatz, reducing the energy of the more obvious 1-parameter square 
structure only slightly. More parameters do not improve the result. The p -condensate lattice is 
superimposed with a more involved lattice of much smaller neutral vortex-like structures, that we 
call superfluid-vortices. 

The superconductivity is mediated by the charged condensate. While we now have proof it is 
organized hexagonally, the structure of the vortices is not essentially different from previous results 
[1,2]. What is interesting is the London relation, which becomes more involved by interaction with 
the neutral condensate. One gets: 

d J 3 -d 3 J = e 2 ml(^ dPl 2 )^ 3 O-l) 

^—do+mk ' 
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Figure 1: Structure of the charged (left) and neutral (right) condensates. 



compared to the result for the ordinary superconductor (2.16) we have that the right hand side 
coefficient is much more 'smoothed out', and in particular does not make the superconductivity 
vanish at the vortex core. 

A second interesting discovery is the existence of the superfluid vortices: configurations with 
a non-trivial winding of the complex phase. The essential difference with the superconducting 
vortices is that the latter tell something about the phase of the charged fields, which is in this case a 
gauge parameter. This allows for the winding property to create a topologically stable lattice. If the 
gauge parameter were not local, we would get a Goldstone boson and see a massless fluctuation, 
rather that seeing it eaten by a photon turned massive. This phenomenon is superfluidity, and can 
be seen in cold liquid Helium by the lack of both translational and rotational resistivity for low 
velocities. 
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Figure 2: Phase of the neutral condensate: left the argument explicitly plotted, right the vorticity of the 
branch points. Red circled are branch points to both p' ' and p. 



The lower order vortices in our superconducting background have indeed a non-trivial winding 
of a non-gauged phase of the neutral fields, hence we call them superfluid. But while the phase 
of p' ' is no symmetry of the DSGS Lagrangian, no massless states are present and there is no 
topological stability. Rather the superfluid vortices are stabilized by the superconducting vortices, 
and we could for example not have an isolated vortex solution. We can see from the plots clearly 
that the superfluid vortices come in pairs, and do not have a net vorticity. 
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4. Outlook 

While we have now computed the structure of the magnetized vacua close to the critical mag- 
netization, an approximation for larger magnetic fields would be interesting. The equations are 
unfortunately much more involved once the condensates get larger. An interesting limit that may 
be in reach could be the case of magnetic field greatly exceeding all mass parameters in the model. 
For the electroweak model this means symmetry breaking is no longer caused by the Higgs field 
and a new phase appears. For the p-mesons no such effect exists, but as the masses become neg- 
ligible, we could approach a phase with more symmetry. The hopes would be then to infer the 
intermediate regime from both near-critical solutions. 

Within the near-critical magnetization considered here, the existence of superfluidity is still 
an unsolved question. No straightforward evidence seems to follow from the theory, but further 
numerical research could shed light on the effect of translations and rotations of the lattice and 
possibly the existence of London-like equations. 

I thank Daniele Binosi for the kind welcome at the TNT-II workshop and Maxim Chernodub 
and Henri Verschelde for collaboration on this subject. 
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